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Abstract 

There is a growing interest in the study of plasma-liquid interaetions with applieation to 
biomedieine, ehemieal disinfeetion, agrieulture, and other fields. This work models the mo¬ 
mentum, heat, and neutral speeies mass transfer between gas and aqueous phases in the eontext 
of a streamer diseharge; the qualitative eonelusions are generally applieable to plasma-liquid 
systems. The problem domain is diseretized using the finite element method. The most in¬ 
teresting and relevant model result for applieation purposes is the steep gradients in reaetive 
speeies at the interfaee. At the eenter of where the reaetive gas stream impinges on the water 
surfaee, the aqueous eoneentrations of OH and ONOOH deerease by roughly 9 and 4 orders of 
magnitude respeetively within 50 qm of the interfaee. Reeognizing the limited penetration of 
reaetive plasma speeies into the aqueous phase is eritieal to diseussions about the therapeutie 
meehanisms for direet plasma treatment of biologieal solutions. Other interesting results from 
this study inelude the presenee of a 10 K temperature drop in the gas boundary layer adjaeent 
to the interfaee that arises from eonveetive eooling and water evaporation. Aeeounting for the 
resulting differenee between gas and liquid bulk temperatures has a signifieant impaet on re- 
aetion kineties; faetor of two ehanges in terminal aqueous speeies eoneentrations like H 2 O 2 , 
NO 2 , and N 03 ^ are observed if the effeet of evaporative eooling is not included. 


1 Introduction 


There is a general interest in the study of plasma-liquid interactions within the scientific community 
for an array of applications, including but not limited to biomedicine and biological disinfection 
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(U|2l|3lm[5113, chemical disinfection ElIHlIll, and agricultural applications. IfTOliTTlI In order 
to sueeessfully realize these applieations and develop mature teehnologies, the basie seience un¬ 
derlying these eoupled plasma gas-liquid systems must be well understood. Recent experimental 
Ell Il3l SH US and modeling work [[T3 [T^ [HI have revealed many of the important physio- 
ehemical processes. However, more work remains, espeeially in the understanding of multi-phase 
eonveetive systems like jets or streamers over water. 

This paper addresses several points regarding transport in plasma-liquid systems. First, eonveetion- 
induced evaporation of water ean lead to rather large temperature gradients at the gas-liquid inter- 
faee, as large as 8 Kelvin/100 microns in the gas boundary layer. The drop in temperature from 
gas to liquid due to evaporative eooling has not yet been considered in the plasma literature; it 
is direetly relevant to plasma ehemistry because of the Arrhenius dependenee of many reaction 
rate constants on temperature. In addition to evaporative eooling, eonveetion leads to depletion of 
water vapor concentration in the vicinity of the diseharge. This will also impact gas-phase plasma 
chemieal reaetions that depend on water vapor as a reactant. |[T9ll 

This study illustrates the sharp drop in the eoncentrations of highly reaetive species from the liquid 
surface to the liquid bulk. In the liquid near the interface, OH eoneentrations fall by as many as 
9 orders of magnitude over tens of mierons. Though the existence of strong gradients in aqueous 
radical concentrations is well known in the bioehemistry literature, Il20l 1^ . it has yet to reeeive 
signifieant attention in the plasma-liquid and plasma-medieine eommunities. Referenee |[T8]I ex¬ 
plicitly models the limited penetration of primary plasma-generated speeies into aqueous solutions. 
Reeognizing these strong gradients in reaetive radical concentrations has important implications 
for understanding plasma therapeutics in which cellular effeets are often observed a signifieant dis¬ 
tance (1 em in ease of subcutaneous tumors Il23 ) away from local plasma treatment. The results 
shown here are eonsistent with the model developed by Graves Il22ll where short-lived plasma¬ 
generated radicals from the gas phase reaet in the surfaee layer of the aqueous/biological phase 
to form longer-lived species such as oxidized and/or nitrated/nitrosylated proteins, peptides, lipids 
whieh ean then enter or oommunieate with surfaee cells. 


2 Model Description 


The eonveetive system ehosen for modelling is shown in fig. It is essentially a point-to-plane 
pulsed-streamer in whieh the liquid surface serves as the gas discharge eathode. The streamer is 
self-pulsed because of a ballasting resistor; typieal diseharge voltages are between 6 and 8 kV 
and pulse frequencies are generally 10 to 30 kHz. The sharp anode tip can be stationed anywhere 
between 3 and 15 mm above the water surfaee. The water is eontained within a glass petri dish of 
radius 3 em; water treatment volumes are generally between 10 and 20 mL. 
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Figure 1: Experimental set-up for pulsed streamer-liquid system. Discharge voltages are typically 
between 6 and 8 kV with pulse frequencies between 10 and 30 kHz. The gap from the needle tip 
to the water surface is between 3 and 15 mm. Water treatment volumes are between 10 and 20 mL. 


We have chosen not to explicitly model the atmospheric discharge due to the computational cost of 
modeling of gas-liquid interfacial transport over time-scales of minutes (the equivalent of millions 
of individual pulses). Instead, gas phase concentrations of the species modeled are taken from 
the DBD-water calculations in [fTTl and diluted by a factor of 10; the gas phase species included 
and their respective inlet concentrations are shown in table The factor 10 dilution is done for 
numerical stability. High concentrations of reactive oxygen and nitrogen species (RONS) at the 
interface combined with very fast rate coefficients, particularly for reactions involving OH, can 
lead to singularities over our time scales of interest unless the computational mesh is sufficiently 
refined, requiring significant computational power. All concentration results and discussion will be 
focused on relative magnitudes and variation with space and time. The results in [fTTl are chosen 
as an input basis because of the absence of experimental work for atmospheric plasmas in which 
a suite of gas-phase plasma species concentrations are reported. Most of the measurements in the 
literature focus on single species like OH [I23ll2^l25ll26l . atomic oxygen [[27l . or NO [l28l that are 
accessible through Laser Induced Fluorescence (LIE). 

A key feature of streamer discharges is the ionic wind. The ionic wind is a net flow of gas towards 
the cathode that results from the effective drag of cations on neutral gas molecules. In the streamer- 
liquid system studied, the ionic wind creates convective forces that affect species transport at the 
gas-liquid interface. Since the discharge is not modelled, another mechanism must be used to 
produce the ionic wind present in a streamer discharge. The mechanism chosen is a jet-like inlet 
with diameter equal to the diameter of the needle used in the experimental streamer set-up. The 
velocity profile at the inlet of the jet channel is chosen such that the channel exit profile closely 
mimics that expected in a streamer experiment. The profile is built from the expected streamer 
maximum axial velocity which is in turn determined from the modeling work of Zhao et. al. Il29ll 
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Using values of 6 kV for the diseharge voltage and 6.5 mm for the gap distanee, the maximum 
axial veloeity is interpolated to be 7.75 m/s near the needle tip. The spatial inlet profile is then 
eomputed from eq. Q assuming laminar no-slip eonditions in the inlet ehannel: [f30l 


U — ^z,niax 




( 1 ) 


The model geometry is shown in fig. model inputs are summarized in tables to The gov¬ 
erning equations used to model the gas and liquid phases are the ineompressible Navier-Stokes 
equations for momentum transport and eonveetion-diffusion equations for heat and mass trans¬ 
port; they are shown in eqs. 0 to Q 


V • M = 0 

(2) 

— +u-Vu j = —V p + iN^u 

(3) 

pCpM-vr = v-(kvr) 

(4) 

7-(-AVQ) + u-VQ = R,- 

(5) 


with u representing the fluid veloeity, p the overall mass density, p the statie pressure, /j the dynamie 
viseosity, Cp the eonstant pressure heat eapaeity, T the temperature, k the thermal eonduetivity, D, 
the diffusivity of speeies i, C, the eoneentration of speeies i, and R,- the souree term representing 
ehemieal reaetions. When solving the model equations, the transient equations]^ andare solved 
until the veloeity eomponents and pressure reaeh a steady-state. Then the steady-state veloeity field 
is inserted into the eonveetion terms in equations and The heat and mass transport equations 
are typieally solved for a physieal time of 1000 seeonds in order to mateh eommon experimental 
treatment times in our laboratory. The temperature dependenee of transport parameters in equa¬ 
tions |^and[^ ineluding reaetion rate eoeffieients, is eonsidered. Gaseous diffusion eoeffieients are 
assumed to seale with T^/^ as predieted by Chapman-Enskog theory. OTl p. 119] The temperature 
dependenee of liquid phase diffusion eoeffieients is eonstrueted using the Stokes-Einstein equation 
[|^ p. 529] and the following equation for viseosity: jSOl p. 31] 


exp{3.STb/T) (6) 

where T^ is the boiling point of the solvent (373 K in the ease of water). Reaetion rate eoeffieients 
and their temperature dependenee are eontained in table Beeause the fluid flow equations are 
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solved prior to solution of the heat transport equation, the temperature dependenee of transport 
parameters in equations and is not ineluded in this work. Model boundary eonditions are as 
follows: for fluid flow, all solid walls are assumed to be no-slip, i.e. all veloeity eomponents at the 
walls are zero. For most results presented the gas-liquid interfaee is assumed to be flat and statie. 
For the statie interfaee, the z-veloeity eomponent at the interfaee is set to zero in both the gas 
and liquid phases. Both the normal and shear stresses are eontinous aeross the interfaee. One may 
question whether the assumption of a flat interfaee is valid when in reality the interfaee is deformed 
by the gas flow. In the results seetion it will be demonstrated that for the purpose of this work, an 
assumption of a flat interfaetion is suffleient. For temperature ealeulations, all solid surfaees are 
assumed to be insulated; at the interfaee, the temperature is assumed eontinuous. Additionally, 
evaporation of water at the surfaee is eoupled to the system’s heat transport in the following way: 
lEQl 


Qb = Jz,H20 ■ Hyap = -DH20,g 


dC, 


H20{g) 


dz 


■H, 


z=interface 


vap 


(7) 


where Q^, is the heat flux, izMiO is the molar flux of H 2 O eoming from evaporation, and Hyap 
is the latent heat of vaporization for water. The eoneentration of water vapor at the interfaee is 
determined from Antoine’s equation: [l32ll 


logioP = 8.07131 — 


1730.63 
233.426-hr 


( 8 ) 


with p in units of mmHg, T in units of degrees Celsius, and the eonstant values taken from [|3^ . 
The temperature at the gas inlet is set to 300 K for all times. Initially the temperature in both gas 
and liquid domains is 300 K. Inlet eoneentrations for speeies other than water vapor are speeified 
at the beginnning of the jet flow development ehannel and are given in table Dilute speeies that 
are present in both gas and aqueous phases (OH, H 2 O 2 , NO, NO 2 , N 2 O 4 , HNO 2 , and HNO 3 ) have 
eontinuous fluxes aeross the gas-liquid interfaee. Coneentrations immediately above and below 
the interfaee are assumed to be in equilibrium as deseribed by their Henry’s law eoeffleients, listed 
in table 1^ A limited set of gas and liquid phase reaetions eneompassing the most important NO;,: 
ehemistry is used to reduee eomputational expense. The reaetions used and their eorresponding 
rate eoeffleients are listed in table It should be noted that beeause of its highly aeidie nature (pKa 
-1.3), HNO 3 is assumed to dissoeiate into H+ and NO 3 immediately after entering the aqueous 
phase. The model equations are solved using the finite element method implemented in Comsol 
Multiphysies version 4.4. Copies of the model are freely available upon request. 

When the “diseharge region" is diseussed, it is in referenee to the gas region roughly demareated 
in the z direetion by the jet outlet/needle anode and the water surfaee, and in the radial direetion 
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by the radius of the jet ehannel/needle, where the plasma is visible during experiments. The reader 
is reminded that the plasma and its eleetrodynamies are not explieitly modeled; the foeus of this 
investigation is the qualitative behavior of momentum, heat, and neutral speeies mass transport in 
eonveetive systems. The qualitative eonelusions drawn here are equally applieable to atmospherie 
jets or streamer systems in whieh the ionie wind plays a key role in momentum transport. More¬ 
over, we postulate that the steep gradients in highly reaetive neutral speeies eoneentrations at the 
gas-liquid interfaee shown in the proeeeding seetion are a universal phenomena of atmospherie 
plasma-liquid systems, regardless of whether they are eonveetive or diffusive; this is eonsistent 
with the reeent researeh by IfTSl as well as the greater reservoir of bioehemistry literature [I2T]| . 
Morever, these gradients have important implieations for plasma medieine, mainly that eellular 
responses must be indueed through seeondary as opposed to direetly generated plasma reaetivity. 


Gas phase speeies 

OH, H 2 O 2 , NO, NO 2 , N 2 O 4 , HNO 2 , HNO 3 , H 2 O 

Liquid phase speeies 

OH, H 2 O 2 , NO, NO 2 , N 2 O 4 , HNO 2 , NOT, NOT, ONOOH, 
H+, OH 


Table 1: Speeies ineluded in the model 


Needle diamater 

1.2 mm 

Petri dish diameter 

6 em 

Gap distanee 

6.5 mm 

Water volume 

10 mL 

Jet ehannel length 

2 em 

Maximum axial veloeity 

7.75 m/s 


Table 2: General model inputs 


3 Results and Discussion 

The steady state veloeity field is shown in figure The reeireulating pattern observed in the 
gas phase is eonsistent with the eorona diseharge modelling results reported in [[29l . A similar 
reeireulation pattern is observed in the liquid phase, indueed by shear stresses between the phases 
at the gas-liquid interfaee. The maximum liquid phase veloeity, . 1 m/s, oeeurs along the interfaee 
approximately .8 mm away from the stagnation point (r = 0). 

The temperature and water vapor profiles at the end of the simulation (t = 1000 seeonds) are 
shown in figures and Notably, the temperature in the bulk liquid has fallen elose to 10 K 
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Molecule 

H, [unitless] 

OH 

6.92-10^ 

H 2 O 2 

1.92-10^ 

NO 

4 . 4 - 10-2 

NO 2 

2 . 8 - 10-1 

N 2 O 4 

3.69-101 

HNO 2 

1.15-103 

HNO 3 

4.8-10^ 

HOONO 

4.8-10^ 


Table 3: The Henry’s constant for a number of molecules ifTTl . 


Molecule 

D [m2 s-i] 

reference 

OH(g) 

4-10-3 

m 

H202(g) 

2-10-3 

m 

NO(g) 

2-10-3 

m 

N02(g) 

1.7-10-3 

m 

N204(g) 

1 -10-3 

m 

HN02(g) 

2.1-10-3 

m 

HN03(g) 

2.1-10-3 

m 

H20(g) 

2.3-10-3 

m 

OH(aq) 

2 - 10 -‘i 

Estimate 

H202(aq) 

1.7-10-‘i 

m 

NO(aq) 

2 . 2 - 10 -‘i 

m 

N02(aq) 

2 - 10 -‘i 

m 

N 204 (aq) 

1.5-10-‘i 

Estimate 

HN02(aq) 

2.5-10-‘i 

By analogy with nitric acid 

HN 03 (aq) 

2.5-10-‘i 

m 

ONOOH(aq) 

2.5-10-‘i 

By analogy with nitric acid 

N02(aq) 

1.7-10-‘i 

m 

N 03 -(aq) 

1.7-10-‘i 

m 

H+(aq) 

7-10-‘i 

m 

OH-(aq) 

1.7-10-‘i 

By anology with other anions 


Table 4: The diffusion coefficients for a number of molecules at 300 K. See text for implementation 
of temperature dependence. 
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Figure 2: Experimental set-up for pulsed streamer-liquid system. Axis units are meters. The 
Python seript used to ereate this figure, as well as the seripts used to ereate all subsequent figures 
ean be found at [I35ll 
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Needle tip/jet outlet 


Petri dish wall 


Gas-liquid interface 
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from its initial value of 300 K. This drop in temperature in the bulk liquid is an example of the 
classical wet-bulb/dry-bulb problem found in chemical engineering texts. Water vapor present in 
the gas above the interface is whisked away by convection. In order to maintain the equilibrium 
vapor pressure required by Antoine’s equation, liquid water must be evaporated, consuming heat 
and lowering the temperature at the interface. Heat then migrates from the bulk liquid to the 
surface, leading to bulk liquid cooling. In this problem we have assumed that the impinging gas 
is at room temperature, and subsequently convection-induced evaporation leads to cooling of the 
bulk liquid below room ambient. However, if there is significant plasma heating of the gas, it is 
possible that heating of the liquid above the initial temperature will be observed experimentally. 
The important point is that the natural coupling between heat transport and evaporation leads to 
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Molecule 

Molecule inlet concentration [m 

OH 

1.3-101*^ 

H 2 O 2 

1 . 6 - 101 ^ 

NO 

00 

0 

00 

NO 2 

5- IQi^ 

N 2 O 4 

0 

HNO 2 

8 - 10 i^ 

HNO 3 

9-10i^ 

H 2 O 

0 


Table 5: Gaseous speeies inlet eoneentrations. IfTTII See text for diseussion 


significant spatial variation in temperature, particularly at the interface, and a cooling of the bulk 
liquid relative to the gas. An example of the steep temperature gradients is shown in figurej^where 
the gas phase temperature changes by 10 K over the span of 200 /jm, immediately above the liquid 
surface. Because reaction rates typically exhibit Arrhenius dependence on temperature, physical 
factors that introduce steep temperature gradients should be included in any model that wishes to 
accurately predict plasma-liquid chemistry. Though not shown here, a simulation was conducted 
in which temperature and water-vapor transport were de-coupled and the interfacial temperature 
gradient removed. Compared to the coupled case, concentrations of long-lived aqueous species 
like H 2 O 2 , NO 2 , and differed by as large as factors of two at the end of the simulation, 
demonstrating the importance of accounting for evaporation-induced temperature gradients. 

Another instance of transport coupling that may impact plasma-liquid chemistry is the significant 
radial gradients in the concentration of water vapor between the needle tip/jet outlet and gas-liquid 
interface. An example of these gradients induced by convective forces are shown in figure 
Examining the dotted curve, which gives the radial distribution of water vapor half-way between 
the jet outlet/needle tip and liquid surface, one can see that the water vapor concentration drops 
precipitously in the region of the discharge (r < 2 mm). In the center of the discharge (r = 0) the 
water vapor concentration is essentially zero. This large drop in the concentration of water vapor 
in the active discharge region due to convection suggests that the concentrations of plasma species 
which rely on H 2 O as a precursor will be reduced at increasing distances from the liquid surface, 
relative to discharges where diffusion is the dominant mechanism of mass transport. 


This model also addresses the role that convection plays in dissolution rates of different gaseous 
species. For this analysis, reactions are turned off, reducing mass transport to only convection and 
diffusion. As one might intuitively expect, the induced convective flow in the liquid significantly 
changes the spatial distribution of aqueous species relative to a diffusion only case, as demonstrated 
for HNO 3 in figures and 8 b However, what is perhaps not intuitive is that though the HNO 3 
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Reaction 

Rate coefficient (Units of s~\ 
moU^ s“\ or mol~^ s^^ Tern- 

Reference 


perature in K. Concentration of M 
in gas and H 2 O in liquid are lumped 
into rate coefficient) 



Gas phase reactions 

1 . 2N02^N204 

2 . N 2 O 4 ^ 2 NO 2 

3. NO + OH + M ^ HNO 2 + M 

4. NO + NO 2 + H 2 O ^ 2 HNO 2 

5. 20H + M^H202 + M 

6 . NO 2 + OH + M -> HNO 3 + M 

7. OH + HNO 2 ^ NO 2 + H 2 O 

8 . 2 HNO 2 ^ NO + NO 2 + H 2 O 

9. HNO 2 + HNO 3 ^ 2 NO 2 + H 2 O 

10. HNO 3 + NO ^ HNO 2 + NO 2 

Liquid phase reactions 

11. N 2 O 4 + H 2 O -> NO 2 + NOg^ + 2H+ 

12. 2 NO 2 + H 2 O ^ NO 2 + NO 3 + 2H+ 

13. NO + NO 2 + H 2 O -> 2 NO 2 + 2H+ 

14. NO 2 + H 2 O 2 + H+ ^ ONOOH + H 2 O 

15. ONOOH ^ .7(N03^ + H+) + .3(N02 + OH) 

16. NO 2 + OH ^ ONOOH 

17. 20H + M^H202+M 

18. HN 02 ^H+ + N 02 


19. H+ + NO 2 -> HNO 2 

20. H+ + OH ^ H 2 O 

21. H20->H+ + 0H 


22. NO + OH ^ HNO 2 

23. 2 HNO 2 ^N0 + N02 


6.02-105.(300/r)3-8 

m 

4.4-10^-exp(-4952/r) 

m 

l.l-10^-(300/r)2-4 

m 

22 

m 

i.o-io^-(r/300)-0'^ 

m 

3.2-10^-(300/r)2-9 

m 

3.0- 106-exp(-390/r) 

m 

6.0-lO-^ 

m 

9.6 

m 

4.4-lO-^ 

m 

1000 

m 

8.4 ■ 10^ ■ exp(-0.033 ■ (T - 295)) 

m 

1.6-105 

m 

1.1-10~5 

m 

.8 

m 

5.3-106 

maizi 

1 -10^ - exp(-450 - (1/r - 1/298)) 

ESI 

3 . 51-106 

From reaction 19 
and equilibrium 
constant 

7-106 

Assumes diffu¬ 
sion limited 

7-106 

Assumes diffu¬ 
sion limited 

7 -10-2 

From reaction 20 
and equilibrium 
constant 

2-10^ 

El 

1.34-10-2.exp(.106-(r-295)) 

m 


Table 6 : Reactions considered in model 
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Figure 3: Velocity magnitude and direction. Axis units are meters. Interface is at z = 0. Axis 
of symmetry is at r = 0. Selected velocity vectors are overlaid on color scale indicating velocity 
magnitude. Red color represents higher values and blue color represents lower values. Velocity 
vector arrows are scaled 12.5 times larger in the aqueous phase. 


spatial distribution changes dramatically depending on whether convection is present in the liquid, 
the volume-averaged uptake of HNO 3 does not change from diffusion-dominated to convection- 
dominated cases as shown in figure If a hydrophobic specie like NO is examined instead of 
a hydrophilic specie like HNO 3 , it is observed that the presence of liquid convection increases 
volume-averaged uptake significantly, as illustrated in figure This fundamental difference in 
behavior between hydrophilic and hydrophobic species can be explained in terms of lumped mass 
transfer resistances. For a hydrophilic specie, the dominant resistance to interfacial transfer is in 
the gas-phase, whereas for a hydrophobic specie the dominant resistance to transfer occurs in the 
liquid phase. Il49l p. 249] Consequently, when convection is added to the liquid phase, effectively 
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r(m) 

Figure 4: 2D temperature profile at t = 1000 seeonds. Red eolor represents higher values; blue 
eolor lower values. Inlet temperature is 300 K. Temperature in the bulk liquid has eooled by 
approximately 10 K beeause of eonveetion-indueed evaporative eooling. 


redueing the liquid-side mass transfer resistanee, the overall resistanee to mass transfer deereases 
and the volume-averaged uptake inereases signifieantly for hydrophobie speeies like NO whereas 
the ehange in overall resistanee is miniseule for hydrophilie speeies like HNO 3 . 

When the full reaetion set is eonsidered, large gradients in reaetive speeies eoneentrations emerge 
at the interfaee. Of partieular interest for biomedieal or pollutant degradation applieations is the 
distribution of hydroxyl radieal in the aqueous phase. Figure [T^ shows a 3D plot of the base 10 log 
of OH(aq) eoneentration versus radial and axial position for t = 1000 seconds. Interfacial gradients 
in OH concentration are prominent in the region where the discharge impinges on the water surface. 
At the stagnation point (r=0) the OH concentration drops by approximately 9 orders of magnitude 
over the span of 50 /jm in the axial direction. Moving away from where the streamer/jet touches the 
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Figure 5: Temperature along z-axis. Illustrates the large temperature gradient that exists at the 
gas-liquid interfaee and the resulting differenee in bulk gas and liquid temperatures. 


surfaee, the interfaeial gradients beeome less pronouneed. Even so, the largest OH eoneentration 
away from the surfaee and in the bulk solution is 3-4 orders of magnitude lower than the peak 
OH eoneentration whieh oeeurs at the interfaee and at the eenter of the impinging streamer/jet. 
The presenee of the liquid phase eonveetive loop is evident from the OH eoneentration hole in the 
eenter of Figure [TO} The effeet is also seen in the eenter of Figure [TT] for ONOOH. 


Another plasma-generated specie which has received much attention in the literature is peroxyni- 
trous acid and its conjugate base peroxynitrite. In the current model, ONOOH exists only in the 
liquid phase and can be created through two mechanisms: reactions 14 and 16 in table [^ Re¬ 
action 14 was the topic of an excellent study in [fT^ and is hypothesized to be a key player in 
the long-term anti-bacterial efficacy of plasma activated water. Figure shows the base 10 log 
of ONOOH concentration as a function of r and z in the aqueous phase. As with OH, there are 
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Figure 6: 2-D water vapor profile at t = 1000 seeonds. Red eolor represents higher values of water 
vapor eoneentration; blue eolor lower values. As implemented through Antoine’s equation, water 
vapor eoneentration at the interfaee is highest where the temperature is highest. Role of eonveetion 
in water vapor profile is evident in deereased eoneentration near the streamer/jet. 


large interfacial concentration gradients. At r=0, the ONOOH concentration drops by 5 orders of 
magnitude over 30 ^um in the axial direction. Away from the stagnation point, especially where the 
convective current flows away from the interface, the interfacial gradient is much smaller but again 
as with OH, the highest bulk concentration is orders of magnitude lower than the peak surface 
concentration. The reason for these large ONOOH gradients is the relative dominance of reaction 
16 over reaction 14 for the given model inputs. Over the course of the simulation, almost 7000 
times more ONOOH is produced through the reaction of OH and NO 2 than through the reaction 
of H+, H 2 O 2 , and NO 2 . As shown in figure [T^ hydroxyl does not penetrate any more than a few 
tens of microns into the liquid phase, so consequently all ONOOH produced through OH is pro- 
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Figure 7: Radial water vapor profiles. Horizontal axis is the radial eoordinate. The top eurve 
(solid) eorresponds to the water vapor eoneentration at the interfaee. The bottom eurve (dashed) 
shows the strong water vapor radial dependenee in the middle of the streamer/jet gap. Gradient is 
largest near and inside the diseharge region. 


dueed within a few tens of mierons of the liquid surfaee. If we examine the produetion of ONOOH 
through reaetion 14, it is observed to be mueh more uniform as illustrated in figure This is 
beeause of the relative uniformity in eoneentration of H+, NO 2 , and H 2 O 2 , although there is a 
peak in their eoneentrations and consequently in their production of ONOOH in the vicinity of the 
impinging streamer/jet. Once the streamer is turned off and surface fluxes of ON and NO 2 are re¬ 
moved, production of ONOOH(aq) will proceed almost exclusively through reaction 14, resulting 
in a mostly homogeneous distribution. This is relevent for applications of PAW since ONOOH is 
the long-term intermediary for production of OH through reaction 15. 

The results presented here suggest two different regimes of activity in the solution: surface and 
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(a) Spatial distribution of HNO3 with liquid convection turned on. 
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(b) Spatial distribution of HNO3 with liquid convection turned off. 
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(a) Comparison of volume averaged uptake of nitrate (HNO 3 before 
dissolution) with liquid convection toggled on or off. Very little 
difference between the two cases. 



(b) Comparison of volume averaged uptake of NO with liquid con¬ 
vection toggled on or off. The presence of convection increases the 
volume-averaged concentration by roughly a factor of 2 over the 
course of the simulation. Note that the vertical scale is much smaller 


for NO than for HNO 3 in 9a because of NO’s hydrophobicity 


Figure 9 
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Figure 10: 3D plot of logio(OH(aq)) as a function of position for t = 1000 seconds. Large interfa¬ 
cial gradients are evident, particularly at the stagnation point (r=z=0). Note the effect of convection 
in the hole in OH concentration in the center of the plot. Some numerical noise is observed at very 
low OH concentrations towards the bottom of the dish. Note that the r and z axes have different 
scales 


bulk. When the discharge is on, reactivity in the form of OH is confined almost entirely to the 
plasma-liquid interface. Bulk solution concentrations of OH are several orders of magnitude lower 
than the surface concentration. This is also true of other potentially biologically important and 
highly reactive reagents such as NO and NO 2 radicals. Less reactive plasma-generated RONS like 
H 2 O 2 and NO 2 have significantly longer lifetimes and can be transported into the bulk solution 
where they can react and form ONOOH, a precursor to OH and NO 2 . However, for the model 
inputs used here, the bulk reactivity generated through H 2 O 2 and NO 2 precursors is orders of 
magnitude less on a per unit time scale than the surface reactivity coming from direct fluxes of 
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Figure 11: 3D plot of logio(ONOOH(aq)) as a function of position for t = 1000 seconds. As with 
OH, large interfacial gradients are evident as is the effect of the liquid convection loop. Note that 
the r and z axes have different scales 


plasma generated radicals. When the discharge is turned off, generation of reactive species through 
H 2 O 2 and NO 2 will persist for a while in the bulk. The results in IIT^ show a nitrite half-life of 2-3 
hours in acidic solution; Traylor et. al. [fTSlI observed anti-bacterial efficacy of PAW, which they 
attributed to ONOOH and its products, for several days after plasma treatment. If one assumes that 
every mole of H 2 O 2 present in solution at the conclusion of our simulation reacts over time to form 
ONOOH and 30% of ONOOH dissociation creates OH, then the amount of OH coming from this 
long-term bulk reaction is 32% less than the amount of OH that comes from surface fluxes while 
the “discharge" is on. This comparison suggests that if one is willing to wait many hours or several 
days, bulk reactivity can approach surface reactivity on an order of magnitude scale. However, 
if an application demands speed and efficient utilization of plasma generated reactivity, the target 
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Figure 12: 3D plot of the base 10 logarithm of the produeation rate of ONOOH through reaetion 14 
listed in tablej^ Beeause of the relative uniformity in the distribution of H+, H 2 O 2 , and NO 2 , the 
produetion of ONOOH through reaetion 14 is mueh more uniform than the overall eoneentration 
profile of ONOOH. Note that the r and z axes have different seales. 


must be plaeed right at the plasma-liquid interfaee. If the target is removed from the plasma by 
millimeters or even hundreds of /jm of aqueous solution, its rate of treatment will rely on a serial 
proeess: the rate of transport of longer-lived speeies like H 2 O 2 and NO 2 through the bulk fluid 
followed by the rate of formation of peroxynitrous aeid chemistry. 

Before concluding, it is worth touching on the assumption of a flat interface. This work discusses 
the role gas phase convection plays in determining momentum, heat, and mass transport across a 
gas-liquid interface and into the bulk liquid. Species uptake is used as the criteria for determining 
whether to include self-consistent deformation of the interface by the gas flow. To make this deter¬ 
mination, the shape of the interface deformation was observed experimentally and the deformation 
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was introduced into the model geometry as shown in figure The depth of the deformation 
was supported by a calculation balancing gravitational (pgh) and convective (^pv^) stresses. The 
fluid flow simulation was then run until it reached stead-state, and then the heat and mass transfer 
equations were solved. In this way, though the interface deformation was not self-consistently de¬ 
termined, its influence on variables of interest could be analyzed. Shown infT^is the effect of the 
interface deformation on total solution uptake of NO(aq). As can be seen in the figure, the effect 
is minute. Because of the characteristics shown in Figure]^ the effect of the interface deformation 
on hydrophilic species transport is expected to be even less. Subsequently, for the purpose of this 
work, the interface deformation can be safely neglected. 
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Figure 13: Geometry for deformed interface simulations. 
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Figure 14: Comparison of total NO(aq) uptake as a function of time for cases in which interface 
deformation is included and not included. As shown in the figure, the interface deformation has 
very little effect on the macrosocpic NO uptake. The effect for hydrophilic species is expected to 
be even less 

4 Conclusion 


A qualitative study of the momentum, heat, and mass transport in a convective plasma-liquid sys¬ 
tem has been conducted. Several interesting results were found. Convective flow of water vapor 
away from the interface leads to a sharp temperature gradient between the bulk gas and liquid 
phases. This sharp gradient is important for accurately determining temperature-dependent rate 
coefficients in both the gas and liquid phases. Additionally, convection drives water vapor away 
from the discharge region except immediately above the liquid surface; this could have important 
consequences for gas-phase generation of reactive species that depend on water as a precursor. 
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Induced convection in the liquid phase substantially changes the spatial distribution of aqueous 
species and increases the volume-averaged uptake of hydrophobic species but interestingly has 
little effect on the volume-averaged uptake of hydrophilic species. This phenomena occurs be¬ 
cause the majority of resistance to interfacial transfer is in the gas phase for hydrophilic species 
and in the liquid phase for hydrophobic species; consequently, decreasing the liquid-phase mass 
transfer resistance by adding liquid convection has a significant impact for hydrophobic but not for 
hydrophilic molecules. 

Perhaps the most interesting result of the study is the sharp distinction between reactivity at the 
surface and in the liquid bulk. Though the limited penetration (tens of /rm or less) of reactive 
neutral radicals is well publicized in prominent biochemistry texts like this phenomena has 
yet to receive significant attention in the plasma-liquid literature with the exception of lITSll . Con¬ 
centrations of species of interest for plasma-medicine and other applications (e.g. OH, NO, NO 2 , 
ONOOH) fall by as many as 9 orders of magnitude within 50 /urn of the interface. In a relatively 
pure aqueous solution as is modeled here, the process responsible for conveying plasma reactivity 
to a target may be transport of H 2 O 2 and acidified nitrite followed by a close-proximity reaction 
to form ONOOH and then OH and NO 2 . For targets in aqueous biological systems, the conveyors 
of reactivity may be proteins modified by primary plasma species. What can be stated from the 
results presented here and in IfTSi is that these conveyors of plasma reactivity are almost probably 
not neutral radicals generated directly by the plasma; rather they are likely secondary species or 
more stable compounds like H 2 O 2 and nitrite. 
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